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Abstract 

Examination of ambient thermodynamic conditions suggest that clathrate hy- 
drates could exist in the martian permafrost, on the surface and in the interior of 
Titan, as well as in other icy satellites. Clathrate hydrates probably formed in a sig- 
nificant fraction of planetesimals in the solar system. Thus, these crystalline solids 
may have been accreted in comets, in the forming giant planets and in their sur- 
rounding satellite systems. In this work, we use a statistical thermodynamic model 
to investigate the composition of clathrate hydrates that may have formed in the 
primordial nebula. In our approach, we consider the formation sequence of the dif- 
ferent ices occurring during the cooling of the nebula, a reasonable idealization of 
the process by which volatiles are trapped in planetesimals. We then determine the 
fractional occupancies of guests in each clathrate hydrate formed at given temper- 
ature. The major ingredient of our model is the description of the guest-clathrate 
hydrate interaction by a spherically averaged Kihara potential with a nominal set 
of parameters, most of which being fitted on experimental equilibrium data. Our 
model allows us to find that Kr, Ar and N2 can be efficiently encaged in clathrate 
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hydrates formed at temperatures higher than ~ 48.5 K in the primitive nebula, in- 
stead of forming pure condensates below 30 K. However, we find at the same time 
that the determination of the relative abundances of guest species incorporated in 
these clathrate hydrates strongly depends on the choice of the parameters of the 
Kihara potential and also on the adopted size of cages. Indeed, testing different 
potential parameters, we have noted that even minor dispersions between the dif- 
ferent existing sets can lead to non-negligible variations in the determination of the 
volatiles trapped in clathrate hydrates formed in the primordial nebula. However, 
these variations are not found to be strong enough to reverse the relative abun- 
dances between the different volatiles in the clathrate hydrates themselves. On the 
other hand, if contraction or expansion of the cages due to temperature variations 
are imposed in our model, the Ar and Kr mole fractions can be modified up to sev- 
eral orders of magnitude in clathrate hydrates. Moreover, mole fractions of other 
molecules such as N2 or CO are also subject to strong changes with the variation 
of the size of the cages. Our results may affect the predictions of the composition 
of the planetesimals formed in the outer solar system. In particular, the volatile 
abundances calculated in the giant planets atmospheres should be altered because 
these quantities are proportional to the mass of accreted and vaporized icy plan- 
etesimals. For similar reasons, the estimates of the volatile budgets accreted by icy 
satellites and comets may also be altered by our calculations. For instance, under 
some conditions, our calculations predict that the abundance of argon in the atmo- 
sphere of Titan should be higher than the value measured by Huygens. Moreover, 
the Ar abundance in comets could be higher than the value predicted by models 
invoking the incorporation of volatiles in the form of clathrate hydrates in these 
bodies. 

1 Introduction 

Clathrate hydrates (hereafter clathrates) were discovered in 1810 by Sir Humphrey 
Davy. Initially considered as laboratory curiosities, it is only from the 1930s that the 
study of their formation conditions became of significant interest because of the clog- 
ging pipelines during transportation of gas under cold conditions. Clathrates are crys- 
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talline solids which look like ice and form when water molecules constitute a cage- 
like structure around small "guest molecules". The most common guest molecules 
in terrestrial clathrates are methane (the most abundant), ethane, propane, butane, ni- 
trogen, carbon dioxide and hydrogen sulfide. Water crystallizes in the cubic system 
in clathrates, rather than in the hexagonal structure of normal ice. Several different 
clathrate structures are known, the two most common ones being named "structure I" 
and "structure II". In structure I, the unit cell is formed of 46 water molecules and can 
incorporate up to 8 guest molecules. In structure II, the unit cell consists of 136 water 
molecules and can incorporate at most 24 guest molecules. 

The thermodynamic conditions prevailing in many bodies of the solar system sug- 
gest that clathrates could also exist in the martian permafrost,— ~2 on the surface and 
in the interior of Titan as well as in other icy satellites. 4 " 8 Moreover, it has been sug- 
gested that the activity observed in some cometary nuclei results from the dissociation 
of these crystalline structures. 9 Generally speaking, clathrates probably participated 
in the formation of planetesimals in the solar system. Indeed, formation scenarios 
of the protoplanetary nebula invoke two main reservoirs of ices that took part in the 
production of icy planetesimals. The first reservoir, located within 30 Astronomical 
Units (AU) of the Sun, contains ices (mostly water ice) originating from the Interstellar 
Medium (ISM) which, due to their proximity to the Sun, were initially vaporized.— 
With time, the decrease of temperature and pressure conditions allowed the water in 
this reservoir to condense at ~ 150 K (for a total gas pressure of ~ 10~ 7 bar) in the 
form of microscopic crystalline ice.— It is then considered that a substantial fraction 
of the volatile species were trapped as clathrates as long as free water ice was available 
within ~ 30 AU in the outer solar nebula. 12 On the other hand, remaining volatiles that 
have not been enclathrated due to the lack of available water ice have probably formed 
pure condensates at lower temperatures in this part of the nebula. 13 ' 14 The other reser- 
voir, located at larger heliocentric distances, is composed of ices originating from ISM 
that did not vaporize when entering into the disk. In this reservoir, water ice was essen- 
tially in the amorphous form and the other volatiles remained trapped in the amorphous 
matrix. 15 ' 16 Consequently, icy planetesimals formed at heliocentric distances below 30 
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AU mainly agglomerated from clathrates while, in contrast, those produced at higher 
heliocentric distances (i.e. in the cold outer part of the solar nebula) are expected to 
be formed from primordial amorphous ice originating from ISM. Thus, clathrates may 
have been accreted in comets, in the forming giant planets and in their surrounding 
satellite systems.*^ ' — 

During the twentieth century, many theoretical and experimental studies allowed 
characterization of the crystalline structures of the most common clathrates. Mean- 
while, a classification has been established to identify the nature of the clathrate and the 
form of occupation of the trapped molecules (single clathrate, multiple guest clathrates, 
etc). From the knowledge of the structure of clathrates, predictive rigorous methods 
have been developed to determine their thermodynamic properties. In particular, van 
der Waals & Platteeuw 25 laid the foundations of a statistical thermodynamics model 
to determine the properties of clathrates. This method is an excellent modern exam- 
ple of the use of statistical thermodynamics to predict macroscopic quantities such as 
temperature and pressure, using the microscopic properties like potential interactions. 
This approach, used today in the industry and in science, has saved substantial experi- 
mental efforts for the determination of i) the equilibrium pressure of a clathrate formed 
from various mixtures and ii) the mole fraction of the different species trapped in the 
clathrate from a given fluid phase. 

In the present work, we use a statistical thermodynamic model derived from the 
approach of van der Waals & Platteeuw in order to investigate the composition of 
clathrates that may have formed in the solar nebula. Indeed, many works published 
in the last decade and detailing the formation conditions of ices in the solar nebula 
have neglected the possibility of multiple guest trapping in clathrates. 14 i 17 i 18 i 20 i 24 i 26 i 27 
In our approach, we consider the formation sequence of the different ices occurring 
during the cooling of the solar nebula and that is usually used to describe the pro- 
cess by which volatiles are trapped in planetesimal s 14 i 17 i 18 i 20 i 24 i 26 i 27 (see Fig. []]). We 
then determine the fractional occupancies of guests in each clathrate formed at given 
temperature. Similarly to papers following the approach of van der Waals & Platteeuw, 
our model is based on the use of intermolecular potentials, which themselves depend on 
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parameters describing the interaction between the molecule and the cage, called "Ki- 
hara parameters". Because the models are extremely sensitive to the choice of these 
parameters, and because different sets of data exist in the literature, 28-34 we examine 
the influence of the interaction potential parameters on our calculations. Moreover, it 
has been shown that the size of the clathrate cages depends on their formation temper- 
ature.—"^ This also leads us to investigate the influence of the size of cages on the 
resulting composition of clathrates formed in the low pressure and temperature condi- 
tions of the nebula. 

The paper is organized as follows. In Section [2] we briefly depict the formation 
sequence of the different ices, including clathrates, in the outer solar nebula. In Section 
[3] we describe the statistical model based on the work of van der Waals and Platteeuw, 
and used to calculate the precise composition of clathrates. In Section|4] we define an 
optimal set of Kihara parameters used in our statistical model in order to calculate the 
nominal fractional occupancies of the different guests incorporated in clathrates formed 
in the solar nebula. Section [5] is devoted to the investigation of the influence of the 
potential parameters and structural characteristics (i.e. size of cages) of clathrates on 
our calculations of the fractional occupancies. In Section|6l we discuss the implications 
of our calculations for the composition of the outer solar system. In particular, we 
focus our attention on the volatile abundances in the atmosphere of Saturn, the issue of 
the measured argon deficiency in Titan and the prediction of the noble gas content in 
comets. Section|7]is devoted to the summary and discussion of our results. 

2 Formation of clathrates in the primordial nebula 

In the present work, we focus our attention on the formation sequence of the dif- 
ferent ices produced in the giant planets formation zone in the primordial nebula. 
Once formed, these ices will add to the composition of the planetesimals accreted 
by the giant planets and their surrounding satellites during their growth. The com- 
position of the initial gas phase of the disk is defined as follows: we assume that the 
abundances of all elements, including oxygen, are protosolar 38 and that O, C, and 
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Table 1 : Gas phase abundances in the solar nebula. 



Species X 


(X/H 2 ) 




species X 


(X/H 2 ) 




O 


1.16 x 10" 


-3 


N 2 


7.62 x 10" 


-5 


C 


5.82 x 10" 


-4 


NH 3 


7.62 x 10" 


-6 


N 


1.60 x 10" 


-4 


CO 


4.91 x 10 


-4 


S 


3.66 x 10" 


-5 


co 2 


7.01 x 10" 


-5 


At 


8.43 x 10" 


-6 


CH 3 OH 


1.40 x 10" 


-5 


Kr 


4.54 x 10" 


-9 


CH 4 


7.01 x 10" 


-6 


Xe 


4.44 x 10" 


10 


H 2 S 


1.83 x 10" 


-5 


H 2 


5.15 x 10" 


-4 









Elemental abundances derive from Lodders (2003). Molecular abundances result from 
the distribution of elements following the approach given in the text. 

N exist only under the form of H 2 0, CO, C0 2 , CH3OH, CH 4 , N 2 , and NH 3 . The 
abundances of CO, C0 2 , CH3OH, CH4, N 2 and NH3 are then determined from the 
adopted CO/C0 2 /CH30H/CH4 and N 2 /NH3 gas phase molecular ratios. Once the 
abundances of these molecules are fixed, the remaining O gives the abundance of 
H 2 0. Concerning the distribution of elements in the main volatile molecules, we set 
CO/C0 2 /CH 3 OH/CH 4 = 70/10/2/1 in the gas phase of the disk, values that are consis- 
tent with the ISM measurements considering the contributions of both gas and solid 
phases in the lines of sight.— "^ In addition, S is assumed to exist in the form of H 2 S, 
with H 2 S/H 2 = 0.5 x (S/H2)gJj, and other refractory sulfide components. 43 We also 
consider N 2 /NH3 = 10/1 in the nebula gas-phase, a value compatible with thermo- 
chemical models of the solar nebula 44 and with observations of cometary comae.— In 
the following, we adopt these mixing ratios as our nominal model of the solar nebula 
gas phase composition (see TableQ}. 

The process by which volatiles are trapped in icy planetesimals, illustrated in Fig. 
Q] is calculated using the stability curves of stochiometric hydrates, clathrates and pure 
condensates, and the thermodynamic path (hereafter cooling curve) detailing the evo- 
lution of temperature and pressure between 5 and 20 AU roughly corresponding to the 
formation locations of the giant planets in the solar nebula |3- We refer the reader to 



'(S/H2) v) means the solar abundance of S relative to tb. is the astronomical symbol for the Sun. 
2 Recent models of giant planets formation show that instead of forming at 20 and 30 AU, i.e. their current 
orbits, Uranus and Neptune underwent most of their growth among proto- Jupiter and proto-Saturn and were 
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Figure 1: Equilibrium curves of NH3-H2O hydrate, H2S, Xe, CH4 and CO clathrates 
(solid lines), CH3OH, CO2, Kr, CO, Ar and N2 pure condensates (dotted lines), and 
thermodynamic path followed by the solar nebula between 5 and 20 AU as a function 
of time, respectively, assuming a full efficiency of clathration. Abundances of various 
elements are solar, with CO/C0 2 /CH 3 OH/CH 4 = 70/10/2/1, H 2 S/H 2 = 0.5 x (S/H 2 )©, 
and N2/NH3 = 10/1 in the gas phase of the disk. Species remain in the gas phase above 
the stability curves. Below, they are trapped as clathrates or simply condense. 
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the works of Papaloizou & Terquem and Alibert et al. for a full description of 
the turbulent model of accretion disk used here. The stability curves of hydrates and 
clathrates derive from Lunine & Stevenson's compilation 4 of published experimental 
work, in which data are available at relatively low temperatures and pressures. On 
the other hand, the stability curves of pure condensates used in our calculations de- 
rive from the compilation of laboratory data given in the CRC Handbook of Chemistry 
and Physics. 50 The cooling curve intercepts the stability curves of the different ices at 
particular temperatures and pressures. For each ice considered, the domain of stability 
is the region located below its corresponding stability curve. The clathration process 
stops when no more crystalline water ice is available to trap the volatile species. Note 
that, in the pressure conditions of the solar nebula, CO2 is the only species that crys- 
tallizes at a higher temperature than its associated clathrate. We then assume that solid 
CO2 is the only existing condensed form of CO2 in this environment. In addition, we 
have considered only the formation of pure ice of CH3OH in our calculations since, to 
the best of our knowledge, no experimental data concerning the stability curve of its 
associated clathrate have been reported in the literature. 

In this study, we assume that the clathration efficiency is total, implying that guest 
molecules had the time to diffuse through porous water-ice solids before their growth 
into planetesimals and their accretion by proto-planets or proto-satellites. This state- 
ment remains plausible only if collisions between planetesimals have exposed essen- 
tially all the ice to the gas over time scales shorter or equal to planetesimals lifetimes 
in the nebula.. 4 . In this case, NH 3 , H 2 S, Xe, CH 4 and - 38.6% of CO form NH 3 -H 2 
hydrate and H2S, Xe, CH4 and CO clathrates with the available water in the outer neb- 
ula. The remaining CO, as well as N2, Kr, and Ar, whose clathration normally occurs 
at lower temperatures, remain in the gas phase until the nebula cools enough to allow 
the formation of pure condensates. Note that, because we assume that the gas phase 
composition of the disk does not vary with the heliocentric distance, the calculated 
clathration conditions remain the same in the 5-20 AU range of the nebula, as shown 

recently by Marboeuf et al.— in the case of their study of the composition of ices pro- 
scattered outward when Jupiter acquired its massive gas envelope, and subsequently evolved toward their 
present orbits.—— 
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duced in protoplanetary disks. Once crystallized, these ices will agglomerate and form 
planetesimals large enough to decouple from the nebular gas and will be accreted by 
the forming planets and satellites. 

3 The statistical-thermodynamic model 

To calculate the relative abundances of guest species incorporated in a clathrate from a 
coexisting gas of specified composition at given temperature and pressure, we follow 
the method described by Lunine & Stevenson 4 and Thomas et al. 3 ' 7 ' 52 which uses 
classical statistical mechanics to relate the macroscopic thermodynamic properties of 
clathrates to the molecular structure and interaction energies. It is based on the original 
ideas of van der Waals & Platteeuw for clathrate formation, which assume that trapping 
of guest molecules into cages corresponds to the three-dimensional generalization of 
ideal localized adsorption. This approach is based on four key assumptions:^! 

1. The host molecules contribution to the free energy is independent of the clathrate 
occupancy. This assumption implies in particular that the guest species do not 
distort the cages. 

2. (a) The cages are singly occupied, (b) Guest molecules rotate freely within the 
cage. 

3. Guest molecules do not interact with each other. 

4. Classical statistics is valid, i.e., quantum effects are negligible. 

In this formalism, the fractional occupancy of a guest molecule K for a given type 
t (t = small or large) of cage can be written as 

Ck,,Pk m 

ykJ i-L,r„/'/- (1) 

where the sum in the denominator includes all the species which are present in the 
initial gas phase. Ckj is the Langmuir constant of species K in the cage of type t , and 
Pk is the partial pressure of species K. This partial pressure is given by Pk = xk x P 
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(we assume that the sample behaves as an ideal gas), with xk the mole fraction of 
species K in the initial gas phase given in Table Q] and P the total gas pressure, which 
is dominated by H2. 

The Langmuir constant depends on the strength of the interaction between each 
guest species and each type of cage, and can be determined by integrating the molecular 
potential within the cavity as 

C,,, = ^/\xp(-^W (2) 
k B T Jo \ k B T ) 

where R c represents the radius of the cavity assumed to be spherical, kg the Boltz- 
mann constant, and wkj (r) is the spherically averaged Kihara potential representing 
the interactions between the guest molecules K and the H2O molecules forming the 

surrounding cage t . This potential w(r) can be written for a spherical guest molecule, 

53 



. W = 2^( 5 -V) + ^"M)-^V) + ^V) 



(3) 



with 



S N (r) = - 



Rc Rc' \ Rc Rc 



(4) 



In Eq. (0, z is the coordination number of the cell. This parameter, which depends 
on the structure of the clathrate (I or II) and on the type of the cage (small or large), is 
given in Table[2] The Kihara parameters a, a and e for the molecule-water interactions 
employed in this work have been taken from Diaz Pefia et al. 30 for CO and from Parrish 
& Prausnitz 28 ' 29 or from Sloan & Koh 34 for all the other molecules of interest (see 
TableEJ. 

Finally, the mole fraction fg of a guest molecule K in a clathrate can be calculated 
with respect to the whole set of species considered in the system as 
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Table 2: Parameters for the cavities. > 



Clathrate structure 



I II 



Cavity type 
Re (A) 



3.975 4.300 3.910 4.730 
2 6 16 8 
20 24 20 28 



small large small large 



b 



z 



R c is the radius of the cavity, b represents the number of small (b s ) or large (b() cages 
per unit cell for a given structure of clathrate (I or II), z is the coordination number in a 
cavity. 



where b s and b\ are the number of small and large cages per unit cell respectively, for 

the clathrate structure under consideration, and with E/V = 1. 

r 

4 Predicted clathrate hydrates occupancies in the pri- 
mordial nebula 

We calculate here the relative abundances of guests that can be incorporated in H2S, 
Xe, CH4 and CO clathrates at the time of their formation in the solar nebula. As far as 
we know, the two most complete sets of Kihara parameters available in the literature 
concerning astrophysical molecules are those published by Parrish & Prausnitz and 
Diaz Pena et al. The parameters given by Parrish & Prausnitz have been obtained by 
comparing calculated chemical potentials based on the structural data of the clathrates 
cages with experimental results based on clathrate dissociation pressure data. These 
parameters have been used in recent work that aimed at investigating the composi- 
tion of clathrates that may exist at the surfaces of Titan and Mars. 2 i 3 i 7 i 52 On the other 
hand, the parameters given by Diaz Pena et al. have been fitted to experimentally mea- 
sured interaction virial coefficients for binary mixtures. These parameters have been 
used to quantify the trapping by clathrates of gases contained in volatiles observed in 
comets. 54 Thomas et al. 52 did a comprehensive comparison between these two sets 



h 



b s yK, s + bjy K ,e 



(5) 
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Table 3: Two different sets for the Kihara potential. 



Set 


Ref. 


Molecule 


G (A) 


e/k B (K) 


a (A) 


(1) 


PP72 


H 2 S 


3.1558 


205.85 


0.36 




PP72 + SK08 


co 2 


2.9681 


169.09 


0.6805 




PP72 


CH 4 


3.2398 


153.17 


0.300 




PP72 


N 2 


3.2199 


127.95 


0.350 




PP72 


Xe 


3.1906 


201.34 


0.280 




PP72 


Ar 


2.9434 


170.50 


0.184 




PP72 


Kr 


2.9739 


198.34 


0.230 




DP82 


CO 


3.101 


134.95 


0.284 


(2) 


SK08 


H 2 S 


3.10000 


212.047 


0.3600 




SK08 


co 2 


2.97638 


175.405 


0.6805 




SK08 


CH 4 


3.14393 


155.593 


0.3834 




SK08 


N 2 


3.13512 


127.426 


0.3526 




SK08 


Xe 


3.32968 


193.708 


0.2357 




PP72 


Ar 


2.9434 


170.50 


0.184 




PP72 


Kr 


2.9739 


198.34 


0.230 




DP82 


CO 


3.101 


134.95 


0.284 



a is the Lennard- Jones diameter, e is the depth of the potential well, and a is the radius 
of the impenetrable core. PP72, DP82 and SK08 correspond to the data taken from 
Parrish & Prausnitz, 28 29 Diaz Pena et al. 30 and Sloan & Koh,— respectively. 

of data and concluded that the parameters given by Parrish & Prausnitz are the most 
reliable because they have been self-consistently determined on experimentally mea- 
sured clathrates properties and also because the results are similar to those obtained 
from other recent compilations (such as the one of Sloan & Koh) whose molecules of 
astrophysical interest are listed in set (2) of Table |3). For these reasons, we have used 
the potential and structural parameters given by Parrish & Prausnitz for H 2 S, CH4, 
N 2 , Xe, Ar and Kr in our nominal calculations. However, in the case of CO, we have 
used the Kirara parameters from Diaz Pena et al. because these data are absent from 
the compilation of Parrish & Prausnitz. Finally, Thomas et al. 3 found that, for com- 
mon molecules, the parameters from Parrish and Prausnitz are very similar to those 
published in recent compilations also obtained from the fit of Langmuir constants to 
simple clathrate formation experimental data,—"— except for the value of the param- 
eter a for C0 2 (the radius of the impenetrable core), which is almost twice larger in 
the recent sets of parameters than in Parrish and Prausnitz's parameters. Following the 
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conclusions of Thomas et al.,^ we have used the CO2 Kirara parameters from Parrish 
& Prausnitz, except for the suspicious aco 2 value, which has been replaced by the one 
given by Sloan & Koh. Set (1) of Table[3]enumerates the list of Kihara parameters for 
various molecules considered in our system and used in our nominal calculations. 

In our calculations, any volatile already trapped or condensed at a higher temper- 
ature than the formation temperature of the clathrate under consideration is excluded 
from the coexisting gas phase composition. This implies that CO2, Xe, CH4, CO, Kr, 
Ar and N2 are considered as possible guests in the case of H2S clathrate. On the other 
hand, only N2, Ar, Kr can become guests in CO clathrate. Figure |2]represents the mole 
fraction / (Eq. 5) of volatiles encaged in structure I and structure II clathrates a priori 
dominated by H2S, Xe, CH4 and CO and formed in the primordial nebula. Interest- 
ingly enough, this figure shows that, contrary to H2S, Xe and CH4, which remain the 
dominating guest species in the clathrates considered, CO becomes a minor compound 
in the clathrate that is expected to be dominated by this molecule. Indeed, whatever 
the structure considered, Kr and N2 become more abundant than CO in this clathrate, 
irrespective of their initial abundances in the gas phase of the nebula. This behavior 
results from the fact that the interaction potential between CO and the water molecules 
forming the surrounding cages is weaker than those involving the other species. 

Table|4]gives the temperature and pressure formation conditions of these clathrates 
in the nebula and the relative abundance /£ defined as the ratio of K/X in X clathrate 
(where X = H2S, Xe, CH4 or CO). Figure [3]represents the abundance ratio F£ which is 
defined as the ratio of K/X in X clathrate (i.e. /£) to K/X in the initial nebula gas phase. 
A guest K which incorporates completely into a given clathrate displays a F£ value of 
1 or greater, provided that there is enough water in the nebula available for clathration. 
This figure shows that volatile species such as Kr, Ar or N2 that are expected to form 
pure condensates below 30 K in the primitive nebula can be efficiently encaged in 
clathrates formed at higher temperature (~ 48.5 K and above). Indeed, irrespective of 
the particular structure, and depending on the amount of water available for clathration, 
these three volatiles can be efficiently trapped in CO clathrate at T = 48.5 K in the 
primitive nebula (panel (d) of Fig. 0). Kr can also be efficiently enclathrated at at T = 
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Figure 2: Mole fraction of volatiles encaged in H2S (a), Xe (b), CH4 (c) and CO (d) 
clathrates. Grey and dark bars correspond to structure I and structure II clathrates, 
respectively. 
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Table 4: Relative abundance f* of volatiles encaged in structure I and structure II 
clathrates. 



Clathrate 




Species 


Structure I 


Structure II 


Jk 


TJ C 




C<J 2 


y .0 / X 1U 


-i 


Q HA w 1 f\- 

8./4X 1U 


-3 


T = 82.3 K 




Xe 


3.28 x 10" 


-6 


7.58 x 10" 


-6 


P = 1.7 x 10" 


-7 u 

bar 


CH4 


8.46 x 10 


-6 


4.51 x 10" 


-5 








9 zL7 v 1 D 
Z.^f / X 1W 


-7 


J.00 X 1U 


-6 






Kr 


2.35 x 10" 


-9 


3.63 x 10" 


-8 






Ar 


4.48 x 10" 


-9 


7.54 x 10" 


-8 






N 2 


4.50 x 10" 


-7 


7.62 x 10" 


-6 


Xe 




CH 4 


1.52 x 10" 


-1 


3.32 x 10" 


-1 


T = 59.8K 
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is defined as the relative abundance of guest K to X in X clathrate (where X = H 2 S, 
Xe, CH4 or CO). Values of T and P correspond to the temperature and pressure of 
the H 2 -dominated gas at which the cooling curve (here at 5 AU) intercepts the equi- 
librium curves of the considered clathrates (see Fig. Q]) Only the species that are not 
yet condensed or trapped prior the epoch of clathrate formation are considered in our 
calculations. 

55.5 K in CH4 clathrate (panel (c) of Fig. [3]). Generally speaking, Figure [3] illustrates 
the fact that some volatiles can be efficiently trapped at relatively high temperature in 
multiple guest clathrates compared to the temperatures at which they are expected to 
condense or to form single guest clathrates in the primordial nebula. It is important 
to note that the present results strongly depend on the choice of the initial gas phase 
conditions and on the amount of water ice available for clathration in the formation 
zone of planetesimals. 
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Figure 3: Abundance ratio F£ is defined as the ratio of K/X in X clathrate (i.e. to 
K/X in the initial nebula gas phase (where X = H 2 S (a), Xe (b), CH 4 (c) and CO (d)). 
Grey and dark bars correspond to structure I and structure II clathrates, respectively. 
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5 Sensitivity to parameters 



The determination of the relative abundances of guest species incorporated in clathrates 
formed in the primitive nebula, and more generally in clathrates formed in any thermo- 
dynamic condition, strongly depends on their structural characteristics (i.e. size of the 
cages) and also on the parameters of the Kihara potential. Here, we investigate the 
influence of these structural characteristics and potential parameters on the fractional 
occupancies of guests in clathrates formed in the nebula. 

5.1 Influence of the interaction potential parameters 

It has been recently shown that, perturbing in the 1-10% range the a and e Kihara pa- 
rameters taken from a given compilation, leads to strong variations of the values of the 
Langmuir constants and thus of the fractional occupancies of enclathrated molecules.— 
In the present case, in order to investigate the sensitivity of the composition of clathrates 
formed in the nebula to the variation of Kihara parameters, we have tested a second 
compilation of parameters close to our nominal set of values taken from Parrish & 
Prausnitz (set 1 of Table but which derives from Sloan & Koh for H2S, CO2, CH4, 
N2 and Xe. Because we did not find any Kihara parameter for Ar, Kr and CO in the re- 
cent published compilations, we have adopted the same values as in our nominal set for 
these molecules (see set 2 of TableO. As a result, the two sets of of Kihara parameters 
used in our comparison are almost identical since several volatiles share the same data 
and the dispersion of data between other compounds is narrow. Two comparisons, rep- 
resented in Figure |U have been made in the cases of Xe and CH4 clathrates produced 
in the primordial nebula. In both cases, calculations have been performed for structure 
I and II clathrates. Despite the similarity of the potential parameters used, one can note 
dispersions up to more than one order of magnitude in the resulting mole fractions (Eq. 
5) for a given volatile encaged in the same clathrate. We then conclude that even minor 
dispersions between the different existing sets of Kihara parameters can lead to large 
variations into the determination of the volatiles trapping in clathrates formed in the 
primordial nebula. However, these variations were not found to be strong enough to 
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modify the relative trapping efficiencies between the different volatiles in the clathrates 
considered here. 

5.2 Influence of cage variations 

In the present work, we have adopted the structural parameters of clathrates given by 
Parrish & Prausnitz and shown in Table [2] Up to now, we have assumed that the size 
of the cages R c is unaffected by the trapping conditions of the different guests in the 
primordial nebula. However, laboratory measurements have shown that the size of the 
cages could increase with temperature and also with the size of the incorporated guest 
species j£ ~ 37 i 56 W e have thus investigated the influence of variations of the cage sizes 
on the mole fractions of guests encaged in clathrates by modifying by up to ±5% the 
values of R c given in Table [2] This large variation is consistent with typical thermal 
expansion or contraction measured in the temperature range 90-270 K.— ~ 37,56 Note 
that because clathration of volatiles occurs at lower temperature in the nebula, slightly 
larger variations of the cage sizes may be expected. 

The evolution of the mole fractions of all guests in structure I clathrates formed 
in the primitive nebula is given in Figure [5] as a function of the size of the cages, 
at clathration temperatures and pressures given in Table [4] This figure shows that 
the contraction or expansion of the cages clearly affects the mole fractions of some 
volatiles in clathrates and that the magnitude of these changes strongly depends on the 
interaction parameters between the guest species and the cages. Indeed, irrespective 
of the clathrate considered, Figure [5] shows that the mole fractions of Ar and Kr can 
vary up to several orders of magnitude by changing the size of the cages in the range 
considered. Moreover, depending on the particular clathrate, mole fractions of other 
molecules such as N2 or CO are also subject to strong changes with the variation of 
the size of the cages. Similar trends have been revealed by performing calculations 
for clathrates of structure II. Our calculations confirm that the trapping propensity of a 
given molecule in clathrates is related rather sensitively to the shape of its intermolec- 
ular potential, which is determined from the adopted Kihara parameters and also from 
the adopted size of the cages (see Eq. [3}. This is to be expected given that the Kihara 
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Figure 4: Mole fraction of volatiles encaged in structure I and structure II clathrates 
dominated by Xe (top panels) and CH4 (bottom panels). Grey and dark bars correspond 
to sets (1) and (2) of parameters of the Kihara potential given in Table [3] respectively. 
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Figure 5: Mole fraction of volatiles encaged in structure I clathrates dominated by H2S 
(a), Xe (b), CH4 (c) and CO (d) as a function of the cage sizes. 
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parameters enter exponentially into the statistical partition function which determines 
the probability of occupancy. The tradeoff between the attractive part of the potential 
and the repulsive part is a function of the guest molecule shape and size, as well as the 
size of the cage, and given that clathrate structures feature two cage sizes, it is not easy 
to intuit what happens to the relative incorporation of a molecule as cages shrink or 
expand thanks to temperature. The statistical mechanical model is therefore essential 
for predictive purposes. 

Figure[5]suggests that thermal variations of the cages need to be taken into account 
in particular if these variations are greater than a few percent. Similar conclusions have 
been obtained by Thomas et al. 52 in their study of clathrate formation and composition 
at the surface-atmosphere interface of Titan. As discussed in the following Section, 
these results may affect the predictions of the composition of the bodies formed in the 
outer solar system. On the other hand, variations with temperature are often not well 
constrained due to the small number of specific systems that have been studied, such 
as for example the pure methane clathrate for which the variations of the cages have 
been found to be small (~ 0.3 %) between 83 and 173 K. 37 In the present study, several 
species are encaged in the same clathrate and the temperature range of interest is lower 
than those considered in experiments. For these reasons, it is difficult to infer variation 
laws describing the expansion/contraction of clathrates at temperatures relevant for the 
solar nebula. 



6 Implications for the composition of the outer solar 
system 

Our calculations of the multiple guest trapping in clathrates formed in the primordial 
nebula have implications for the formation and composition models of the giant planets, 
their surrounding satellites and also comets. Indeed, the fact that several compounds 
expected to be trapped or condensed at low temperature are incorporated at relatively 
higher temperature in clathrates in the nebula modifies the predictions of the composi- 
tion of the icy planetesimals from which these bodies were presumably formed relative 
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to what is obtained from simple condensation models. A few examples are discussed 
here. 

6.1 Volatile enrichments in Saturn 

Measurements by the mass spectrometer aboard the Galileo probe have shown that 
the abundances of C, N, S, Ar, Kr and Xe are all enriched by similar amounts with 
respect to their solar abundances in the atmosphere of Jupiter. 15 i 57 i 58 Similarly, re- 
cent Cassini CIRS observations have also confirmed what was inferred from previous 
measurements, that C is substantially enriched in the atmosphere of Saturn. 59 ' 60 In or- 
der to interpret these volatile enrichments, it has been proposed that the main volatile 
compounds initially existing in the solar nebula gas phase were essentially trapped by 
crystalline water ice in the form of clathrates or hydrates in the feeding zones of Jupiter 
and Saturn. 14 i i7 r 1 9 1 24 1 26 1 27 These j ces then agglomerated and formed planetesimals that 
were ultimately accreted by the forming Jupiter and Saturn. This is then the fraction of 
these icy planetesimals that vaporized when entering the envelopes of the two growing 
planets which engendered the observed volatile enrichments. 

On the other hand, our statistical model allows us to infer that multiple guest 
clathrates were more likely formed than single clathrates in the solar nebula, imply- 
ing substantial changes in the presumed composition of the planetesimals formed from 
these ices. As a result, the volatile enrichments calculated in the giant planets atmo- 
spheres should also be altered because these quantities are proportional to the mass of 
accreted and vaporized icy planetesimals. Indeed, for instance, a recent interpretation 
of the carbon abundance in the atmosphere of Saturn is based on the hypothesis that the 
giant planet never formed at a disk temperature below 30 K, 27 implying that the plan- 
etesimals accreted by the giant planet were impoverished in Ar, Kr, CO and N2 (i.e. 
volatiles whose condensation curves are located below 30 K in the nebula - see Fig. 
Q}. However, our nominal calculations predict that Kr can be entirely trapped in CH4- 
dominated clathrates at ~ 55 K in the nebula, provided that there is enough water avail- 
able for clathration at this temperature in the feeding zone of Saturn. In this scenario, 
the trapping of Kr at higher temperature in the planetesimals accreted by Saturn implies 
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that the atmospheric abundance of this noble gas should also be enhanced by an amount 
similar to that of carbon compared to solar (at least ~ 9 times solar 14 - 60 ), instead of 
being predicted in solar abundance in Saturn's envelope. 27 However, this conclusion is 
valid only for the structural characteristics and the Kihara parameters adopted in our 
nominal calculations. As shown in Figure [3] any expansion of the cage sizes during 
the formation of CPU-dominated clathrate in the solar nebula could strongly reduce the 
trapping efficiency of Kr in this clathrate and thus decrease its resulting enhancement 
in the atmosphere of Saturn. 

6.2 The argon deficiency in Titan 

A puzzling feature of the atmosphere of Titan is that no primordial noble gases other 
than argon were detected by the Gas Chromatograph Mass Spectrometer (GCMS) 
aboard the Huygens probe during its descent to Titan's surface on January 14, 2005. 
The observed argon includes primordial 36 Ar, i.e. the main isotope, and the radiogenic 
isotope 40 Ar, which is a decay product of 40 K— In any case, the 36 Ar/ 14 N is lower 
than the solar value by more than five orders of magnitude. 61 The other primordial 
noble gases Kr and Xe (and 38 Ar) were not detected by the GCMS instrument down to 
upper limits of 10 parts per billion relative to nitrogen. 61 The Kr and Xe deficiencies 
could simply be explained by the presence of clathrates on the surface of Titan that 
would have efficiently incorporated these noble gases. 7 - 52 On the other hand, in order 
to interpret the Ar deficiency in the atmosphere of Titan, it has been proposed that the 
satellite was formed from icy planetesimals initially produced in the solar nebula and 
that were partially devolatilized at a temperature not exceeding ~ 50 K during their 
migration within Saturn's subnebula. 13 In this case, because Ar is poorly trapped in 
clathrates formed above ~ 50 K in the nebula, only tiny amounts of this compound 
would have been incorporated in the building blocks of the forming Titan, in agree- 
ment with the observations. In particular, our nominal model predicts that this noble 
gas would remain essentially trapped in CH4-dominated clathratq 3 ], and subsequently 
in the satellite, with A1/CH4 of 2.7 x 10~ 5 (see Fig. |3J, in good match with the abun- 



The amount of Ar found in Xe and FbS-dominated clathrates is negligible. 
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dances of Ar observed in Titan's atmosphere. However, this statement remains valid 
only if one considers our nominal calculations of multiple guest trapping in clathrates 
formed in the feeding zone of Saturn. Indeed, Figure [5] shows that the trapping of Ar 
can strongly increase by a factor of more than 300 in CH4-dominated clathrate formed 
in the solar nebula if the size of the cages decreases by a few percent. In this case, 
the amount of argon predicted by our calculations in the atmosphere of Titan would be 
higher than the value measured by Huygens or an alternative scenario must be invoked 
to explain its apparent depletion. 

6.3 Noble gas content in comets 

It has been proposed that the composition of volatiles observed in comets could be 
explained on the basis of their trapping in the form of clathrates in the primordial 
nebula. 54 In this model, the key parameter to explain the volatile content observed in 
comets is the amount of water ice available in the region of the nebula where the clathra- 
tion took place. The mass of available water is then varied between the amount needed 
to trap the whole mass of volatiles present in the nebula (high-mass ice scenario) and 
that which is just required to enclathrate volatiles at temperatures above ~ 50 K in the 
nebula (low-mass ice scenario). In this model, the case of argon is interesting because 
its relative abundance in comets is found to be very small (Ar/H20 ~ 10~ 8 ) when the 
low-mass ice scenario is considered. On the other hand, if enough water ice is present, 
Ar/H20 jumps to potentially detectable values (AJ/H2O ~ 10~ 4 -10~ 3 ). However, as 
in the case of Titan described above, taking into account the uncertainties of our cal- 
culations (in particular on the cage sizes), substantial amounts of Ar can be trapped in 
CH4-dominated clathrates. In this case, we find A1-/H2O - 10~ 5 -10- 4 in the low -mass 
ice scenario, a value which is close to the one determined in the case of the high-mass 
ice scenario. Our estimate for the low-mass ice scenario implies then that it is diffi- 
cult to consider the measurement of the abundance of argon in comets as a key test to 
constrain the mass of water ice that was available in the nebula for forming comets. 
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7 Summary and discussion 

In this paper, we have calculated the relative abundances of guest species that can be 
incorporated in clathrates formed in the gas phase and thermodynamic conditions of 
the primordial nebula. We have assumed that the clathration efficiency is total in the 
primitive nebula, implying that guest molecules had the time to diffuse through porous 
water-ice solids before their growth into planetesimals and their accretion by proto- 
planets or proto-satellites. This statement remains plausible only if collisions between 
planetesimals have exposed essentially all the ice to the gas over time scales shorter or 
equal to planetesimals lifetimes in the nebula. 4 However, it is important to note that the 
efficiency of collisions between planetesimals to expose all the "fresh" ice over such a 
time scale still remains questionable and that we have no evidence that clathration was 
important in the primordial nebula. 

The results presented here derive from the usual statistical model based on the van 
der Waals and Platteeuw approach, generalized by Parrish & Prausnitz for the calcu- 
lations of dissociation pressures of multiple guest clathrates. The major ingredient of 
our model is the description of the guest-clathrate interaction by a spherically averaged 
Kihara potential with a nominal set of parameters, most of which being fitted on ex- 
perimental equilibrium data. Our model allows us to find that Kr, Ar and N2 can be 
efficiently encaged in clathrates formed at temperatures higher than ~ 48.5 K, the tem- 
perature of CO clathrate, in the primitive nebula, instead of forming pure condensates 
below 30 K. However, we find at the same time that the determination of the relative 
abundances of guest species incorporated in these clathrates strongly depends on the 
choice of the parameters of the Kihara potential and also on their adopted structural 
characteristics. Indeed, testing different potential parameters, we have noted that even 
minor dispersions between the different existing sets can lead to non negligible varia- 
tions in the determination of the volatiles trapped in clathrates formed in the primordial 
nebula. However, these variations are not found to be strong enough to reverse the rela- 
tive abundances between the different volatiles in the considered clathrates. Moreover, 
we have found that the contraction or expansion of the cages, due to temperature vari- 
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ations, can alter the mole fractions of some volatile molecules up to several orders of 
magnitude in clathrates. On the other hand, due to the lack of laboratory experiments 
describing properly the variation of the size of cages as a function of temperature for 
molecules of astrophysical interest, it is actually difficult to quantify the influence of 
this parameter on the composition of planetesimals formed in the nebula. The deter- 
mination of specific laws describing the variation with temperature of the size of cages 
is important because it could help to constrain the thermodynamic conditions encoun- 
tered by planetesimals during their formation in the primordial nebula. Indeed, given 
a sufficient number of species whose abundances are determined in a particular object, 
our statistical mechanical model could be used to predict the composition and temper- 
ature (through the thermal expansion/contraction of the ice which affects cage size) of 
the planetesimals from which the object formed and hence the composition of bodies 
issued from the same family. Such a measurement will be achievable by the Rosetta 
spacecraft mission towards Comet 67P/Churyumov-Gerasimenko or by the next gen- 
eration of missions designed to investigate the outer solar system. 

Our results alter the predictions of the composition of the planetesimals formed in 
the outer solar system. In particular, the volatile abundances calculated in the giant 
planets atmospheres should be recalculated because these quantities are proportional 
to the mass of accreted and vaporized icy planetesimals. For the same reasons, the esti- 
mates of the volatile budgets accreted by icy satellites and comets must be reconsidered 
based on our calculations. For instance, our calculations predict that the abundance 
of argon in the atmosphere of Titan might for some initial conditions be higher than 
the value measured by Huygens because substantial amounts of this volatile may be 
trapped in the satellite's building blocks during their formation. Similarly to the case 
of Titan, if comets were agglomerated from clathrates in the nebula, then the Ar abun- 
dance in these bodies should still remain potentially detectable, irrespective of their 
formation temperature (~ 50 K and below) because substantial amounts of this volatile 
are expected to be trapped in CH4-dominated clathrate at ~ 55.5 K. 

It is important to note that the validity of the statistical model to determine the 
composition of clathrates in the nebula probably meets some limitations. In particular, 
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the transferability of the Kihara parameters to temperatures and pressures beyond the 
range for which they have been fitted is uncertain. 62 ' 63 Moreover, inconsistencies have 
been evidenced between Kihara parameters derived from different sets of experimental 
data.— In these conditions, supplementary measurements of the equilibrium pressure 
of clathrates at low temperature would be needed in order to obtain a set of Kihara 
parameters consistent with the low temperature and pressure conditions of the nebula. 
Moreover, it has been claimed in the literature that Kihara potential may not accurately 
describe the interaction between guest and water molecules. 65 Recent works have thus 
been based on an atom-atom description of the intermolecular guest-clathrate interac- 
tions, in which effective parameters for these interactions have been fitted from results 
of ab initio quantum mechanical methods.—"— However, the accuracy of this atom- 
atom approach is strongly dependent on the ingredients of the fit, which can be the 
number of water molecules considered in the ab initio calculations, the number of sites 
chosen to represent the interacting molecules, the level of accuracy of the quantum 
methods used... It is then unfortunately very difficult to discriminate which approach 
(atom-atom or spherically averaged) is the most accurate. 
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